RATIONAL APPROXIMATION 
FOR A QUASILINEAR PARABOLIC EQUATION 

P. M. GAUTHIER AND N. TARKHANOV 

Abstract. Approximation theorems, analogous to known results for 
linear elliptic equations, are obtained for solutions of the heat equation. 
Via the Cole-Hopf transformation, this gives rise to approximation the- 
orems for a nonlinear parabolic equation. Burgers' equation. 
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1. Introduction 

It is well known that each holomorphic function on a compact set K in 
the complex plane can be approximated by rational fractions with poles 
away from K. In fact the set of poles can be fixed arbitrarily to satisfy the 
only condition that it meets each connected component of the complement 
of K. If it has at least one limit point at each component of C \ then 
it is actually possible to approximate by finite linear combinations of the 
Cauchy kernel. 

The Cauchy kernel is a fundamental solution for the Cauchy-Riemann 
operator in C. This gives insight into what can be though of as rational 
solutions to linear partial differential equations. These are just finite linear 
combinations of a right fundamental solution to the given equation. In 
this way the rational approximation theory naturally extends to solutions 
of linear elliptic equations, see [lU 5.3.2]. 
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Since the analogue of rational approximation has been well developed for 
elliptic equations, it is natural, from the theoretical viewpoint, to attempt 
a similar theory for parabolic equations. Thus, we seek to approximate a 
solution of a parabolic equation by solutions having isolated singularities, 
and we also attempt to specify the location of these singularities. 

As a possible physical interpretation for rational parabolic approximation, 
consider the following. Let B be an open set in M" and let n be a heat 
distribution on B, that is a solution of the heat equation on M x S C M""''^. 
We call the pair (B, u) a thermal box. Thus, a thermal box is a domain 
endowed with a heat distribution. We may think of ovens and refrigerators 
as examples, but we may also think of a house or a region of physical space 
as examples. In the case that S is a refrigerator, an oven or a house, we are 
usually interested in prescribing a temperature distribution u{t,x), which is 
constant as a function of x. But we could ask for much more. We could ask 
for different temperatures in different parts of B and, indeed, this is often 
done in houses. For example, let Bi and i?2 be two subregions of B and let 
ui and U2 be heat distributions on Bi and B2, respectively. Is it possible to 
design the thermal box (B, u) to these specifications? That is, can we find a 
temperature distribution u on B such that u = ui on Bi and u = U2 on B2"! 
Since the heat distribution u is analytic in the space variable, the answer is in 
general "no." However, from the point of view of every conceivable physical 
application, it would be sufficient to find on B, which approximate ui 
and U2 on Bi and B2, respectively, within some tolerance e, provided we 
can do this for arbitrarily small tolerances e. This is precisely the kind of 
approximation we wish to investigate: uniform approximation on subsets of 
B by solutions of the heat equation on all of B. From the engineering point 
of view, one is interested not merely in the mathematical existence of the 
global approximating heat distribution u^. One seeks a way of designing 
B to obtain such a heat distribution u^. A mathematical model for this 
is to take as u a finite linear combination of fundamental solutions whose 
singularities lie on the boundary of B. We shall call these simple rational 
functions. We may think of such a boundary point as a heat source, or 
infinitely hot spot, if the coefficient is positive at that point, and we may 
think of it as a heat sink or infinitely cold spot, if the coefficient is negative. 
As an engineering approximation, one can impose hot temperatures at the 
hot spots and cold temperatures at the cold spots on the walls of the thermal 
box B. 

In addition to proving results on approximation by rational solutions to 
parabolic equations, we shall also investigate the problem of approximating 
solutions in a one domain by solutions in a given larger domain. Pairs of 
domains where this is possible are called Runge pairs. 

In 1929, Birkhoff [2] showed the existence of an entire function u{z) with 
the following remarkable property. For any entire function f{z), there is a 
sequence {aj} such that the sequence {u{z+aj)} of translates of u converges 
to the function f{z) uniformly on compact subsets of C. Thus, the translates 
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of the function u approximate all entire functions /. Such a function u is 
said to be universal. 

Another kind of universality is that of universal series. Let V he a, space 
of functions and {uj} a sequence of functions in V. The series Yl'^^j is 
said to be universal in V if the set formed by the partial sums is dense 
in V. Of course, such series are to be thought of as formal series. They 
are never convergent, if V is not trivial. In 1951, Seleznev [15] showed 
the existence of a universal power series and this result was refined in ^12j . 
[3] and [13]. Subsequently, other universal series were found for solutions 
to the Cauchy-Riemann as well as the Laplace equation and other elliptic 
equations. Recently, an abstract theory of universal series, covering most of 
these results, were introduced in [T3]. For a broader view of universality, we 
refer the reader to the excellent survey [6]. 

This work is intended as an attempt to specify these results in the con- 
text of nonlinear partial differential equations. We shall present results on 
universal functions as well as universal (formal) series for solutions to para- 
bolic equations. The crucial role in our investigation is due to the so-called 
Cole-Hopf transformation. 



The Cole-Hopf transformation was discovered independently by Cole [1] 
and Hopf [8] around 1950. It changes Burgers' equation ut + uUx = Uxx 
into the heat equation vt = Vxx- To derive the transform, we let u = —px- 
Then Burgers' equation can be integrated yielding pt — = pxx up to 
a function depending on t only. Let p = 21ogf. Applying some algebra to 
this we get vt = Vxx- 

The n -dimensional forced Burgers equation ut+{u-'S/)u = i^V^u— V/(i, x) 
for u = — Vp, which describes the dynamics of a stirred, pressure less and 
vorticity-free fluid, has found interesting applications in a wide range of 
non-equilibrium statistical physics problems, see [1]. Here, v stands for the 
viscosity. The associated Hamilton-Jacobi equation, satisfied by the velocity 
potential p, 



has been frequently studied as a nonlinear model for the motion of an in- 
terface under deposition, when the forcing potential / is random, delta- 
correlated in both space and time. Equation ()2.ip is the well-known Kardar- 
Parisi-Zhang equation, see pLOJ. 

Starting with this example, we consider a quasilinear partial differential 
equation 



in M""*" , where o is a continuous real- valued function on the real axis. 
Choose a strictly monotone increasing function u = U{p) on M with 



2. The Cole-Hopf transformation 






Pt = Ap + a(p)|Vp 
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the property that 
for all p eR. 

The general solution of this ordinary differential equation satisfying the 
initial condition U'{0) = C/i > is 

(2.3) U'{p) = C/i exp - y a{t)dt^ , 

which is a smooth function on M with positive values. The function u = U{p) 
may be found by integration. In this way we recover what is referred to as 
the Cole-Hopf transformation. 

A simple computation shows that the change of variables u = U{p) re- 
duces (j2.2p to the heat equation 



(2.4) ut = Au 

for the new unknown function u. Hence, the general solution to (j2.2p is 
p{x) = U~^{u{x)), with u satisfying ([27 



Example 2.1. Let a be constant. Then 

\ TT rr 1 ~ exp(— ap) 
a 

U-\u) = -ilogfi 



a ° V Ui 

Using the function U allows one to endow the set of solutions to equation 
(j2.2p with the symmetry pi op2 := U~^{U{pi) + U{p2))- 

3. Rational approximation 

If U is an open subset of we shall denote by S{U) the family of 

all complex- valued solutions u G C°°{U) of the heat equation ut = Au on 
U. The topology on S{U) induced by embedding this space into C°°{U) 
is actually equivalent to the topology of uniform convergence on compact 
subsets of U. 

If S is an arbitrary subset of R"+-^, we shall denote by S{S) the family of 
germs on E of solutions u to the heat equations on some open set (depending 
on u) containing E. Such solutions form a complete locally convex space 
under the inductive limit topology. 

Functions in 5(]R"+^) will be called entire solutions of the heat equation. 
They have singularities at points at infinity of a suitable compactification of 

Set 



— 77- exp — , if t > 0, 

(47rt)"/2 ^ V 4t / ' 

0, if t < 0. 
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This function is locally integrable over M"^ , infinitely differentiable in 
W"'^^ \ {0}, and it satisfies {dt — A)<1> = 6 in the sense of distributions 
in M*^^ , 6 being Dirac's measure at G W'-+'^. The distribution <P IS re- 
ferred to as the standard fundamental solution of convolution type for the 
heat equation on 

Solutions of the heat equation of the form 
J 

(3-1) E E ^3,o.dld'^<l>{t -to,x- xo), 

j=0 \a\<A 

where cj^a G play the role of rational solutions with pole at the point 
{to,xo) £ M""*"^. These are nothing but the potentials of distributions sup- 
ported at (to, Xo). 

By a simple rational solution of the heat equation on R"-"*"^ we mean any 
finite linear combination 

N 

y^C^ 'Pit -tu,X - Xu) 
v=\ 

of the heat kernel itself with poles at points (ti, xi), . . . , (tAr, xat). From this 
point of view, the following is an analogue, for the heat operator, of Runge's 
theorem on approximation by partial fractions for the Cauchy-Riemann op- 
erator. 

Theorem 3.1. For each compact set K C M""*"^, the simple rational solu- 
tions to the heat equation with poles outside of K are dense in the space 
S{K). 

Proof. Let u € S{K). We may suppose that u G C^iW"-^^) has compact 
support and satisfies dtu = Au in a neighbourhood U of K. Thus, 

= /.(*-t'.x-xO(a,„-^„)(*',.Vt'<^' 

= / ^{t- t', X - x') {dtu - Au) {t', x')dt'dx', 

J supp (dtu— Au) 

where dt'dx' denotes Lebesgue measure on M""*"^. For fixed {t,x) away from 
the support of dtu — Au, this is a Riemann integral. 

In order to estimate this integral by Riemann sums, let Q be a compact 
set containing the support of dtu — Au and disjoint from K, such that Q 
is a finite union of hypercubes whose sides are parallel to the coordinate 
hyperplanes. The integral over supp (dtu — Au) is the same as the integral 
over Q, and it is easier to take Riemann sums over Q. liV = {Qi, . . . ,Qn} 
is a partition of Q into hypercubes, we denote by \V\ the mesh of V, that is, 
the maximum of the diameters of the hypercubes of the partition V. For each 
Qi, of the partition V, choose a point {tu,Xy) G Qy. For each (t,x) £U\Q, 



Serge Lange, near the end of his life, asserted in a lecture at the Universite de Montreal, 
that the heat kernel is the most important object in all of mathematics. 
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denote by S-p{t,x) the corresponding Riemann sum of the integral over Q. 
Then, for any fixed {t, x) G U\Q, the Riemann sums converge to u{t, x), as 
\V\ 0. 

Each such Riemann sum has the form 

N 

S-pit, x) ='^c„ <P{t -t„,x - x^), 

where the points (tu,Xi,) are in Q. This Riemann sum IJ-p{t,x) is a simple 
rational solution to the heat equation. Thus, we have simple rational solu- 
tions to the heat equations Up which converge pointwise to u on U \Q, as 
\V\ 0. Since 

{t, x; t', x') ^ <P{t -t',x- x') {dtu - Au) (t', x') 

is uniformly continuous of (U \ Q) x Q, the convergence of this Riemann 
sums (rational solutions) to u is in fact uniform on U \ Q, as {Vl 0. 

It is not hard to see that not only do these rational solutions to the heat 
equation converge to u, but also their partial derivatives converge uniformly 
on U \ Q to the corresponding partial derivatives of u. Thus, we in fact have 
C°° approximation. □ 

One sees that the proof actually goes through for solution of any differen- 
tial equation possessing a left fundamental solution smooth away from the 
diagonal. 

The following result can be thought of as a theorem on rational approx- 
imation with a fixed set of singularities. It contains Theorem 13.11 as a very 
particular case, however, the proof is no longer constructive. 

Theorem 3.2. Let K be a compact set in IR"'+^ and p be a subset ofW^~^^\K 
with the property that any solution g to —dtg = Ag in W^^^ \ K which, 
together with its derivatives in t up to order J and in x up to order A, 
vanishes on p, is zero in some layer around K . Then the potentials 1^3. 
where (to , 2:0) £ P O'^d cj^a £ cli^g dense in the space of all solutions to the 
heat equation on K. 

By a layer around K is meant any open set U \ K, where C/ is a neigh- 
bourhood of K in 

Proof. Denote by S{K) the space of all solutions of the heat equation on 
K. Since the heat equation is hypoelliptic, each distribution u satisfying 
(|2a on an open set U C M"+^ belongs actually to C°°{U). Hence S{K) 
can be specified as a closed subspace of C^{K), this latter is the space of 
all functions on neighbourhoods of K endowed with the inductive limit 
topology. 

Let TZ stand for the subspace of S{K) consisting of the potentials (j3.ip . 
where {tQ,xo) G p and Cj a G Our task is to show that TZ is dense in 
SiK). 
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To this end we use the Hahn-Banach theorem. Pick a continuous hnear 
functional T on C°°(Er). Since the dual space for C°°(i^) just amounts 
to E'^^ the space of all distributions on M""*"^ with support in there is a 
V G E'j^ such that T(u) = (w, u) for all u G C°°(i^). Consider the convolution 
of distributions 

g{i!,x') = {^{t-t',x-x'),v{t,x)) 
= (p' *v, 

where 'P'{t,x) = ^{—t, —x) is the transposed kernel. 

Since ^' is the fundamental solution of the transposed heat equations on 
M"+\ it follows that 

{-dt'-A^,)g = {-dt> - A^,)^' *v 
= 6 * V 
= V 

on all of M""*"^. In particular, g' is a solution of the transposed heat equation 
in the complement of K. 
If J- vanishes on TZ, then 

di,d^,g{t',x') = {di,d^,^{t-t',x-x'),v{t,x)) 

= i-iy+\'^\:F[did^'P{--t',--x')) 
= 

for all {t',x') € p and for all j and a satisfying j < J and |a| < A. By 
assumption, there is a neighbourhood U of -fC in M""'"^, such that g = in 
U\K. 

Thus, if M is a solution of the heat equation in a neighbourhood of K, 
then we obtain 

T{u) = {{-dt-A)g,u) 
= {g,{dt-A)u) 
= 0. 

Now the assertion follows from the Hahn-Banach theorem. □ 

4. Choice of pole sets 

To effectively use Theorem 13.21 we should be able to show explicit sets 
p for which the hypotheses of the theorem are fulfilled. In appropriate 
sense p G R'"+^ \ K are uniqueness sets for solutions of the transposed 
heat equation in the complement of K. More precisely, each solution g to 
~dtg = Ag in M"+^ \ K satisfying dfd^g = on p for all j and a with 
j < J and \a\ < A must vanish in some layer around K. Such results 
can be derived from the analyticity of solutions to the heat equation on 
characteristics. 
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We first recall an idea of [llj wfiicli is known as elliptic continuation of 
solutions of a parabolic equation. 

Theorem 4.1. Let g he a hounded solution to dtg + Ag = in a cylinder 
Z = (0,to] ^ B{0,R). Then for each R' < R there are positive constants e 
and C depending on to, R and R' , such that in Z' = (to — e, to + e) x -6(0, R') 
there is a solution g' to the elliptic equation dfg' + Ag' = with Cauchy 
data 

g'{to,x) = g{to,x), 
dtg' {to, x) = 

for X G -6(0, R'), satisfying 

sup \g'{t,x)\<C sup \g{t,x)\. 
{t,x)eZ' {t,x)ez 

Proof. The existence follows from the Cauchy-Kovalevskaya theorem once 
we observe that the restriction g to the characteristic hyperplane t = to is 
real analytic in x G B{0,R'). The estimate is a consequence of an explicit 
construction. □ 

From Theorem 14. II we deduce in particular that if (7 is a bounded solution 
to dtg + Ag = in the cylinder Z and, for some xo G -6(0, R), the restriction 
g{to,x) decreases faster than any power |x— xo|*^ as x ^ xq, then g{to,x) = 
for all X E B{0,R). 

Example 4.2. Set p = U\K where U is an arbitrary neighbourhood of K 
in M""*"^. Then the hypotheses of Theorem 13.21 are obviously fulfilled with 
J = and A = 0. 

Example 4.3. Let U be an arbitrary neighbourhood of K in M""'"^ and p 
a subset oi U \ K, such that, for any hyperplane H orthogonal to the time 
axis, each connected component oi H f] {U \ K) contains at least one point 
of p. By the above, the hypotheses of Theorem 13.21 are fulfilled with J = 
and A = 00. 

Example 4.4. Suppose K is a compact set, such that each hyperplane 
orthogonal to the time axis meets K in a set with connected complement. 
Let p be any curve in M""'"^ \ K whose projection to the time axis contains 
that of a neighbourhood of K. Then the hypotheses of Theorem 13.21 are 
fulfilled with J = and ^4 = 00. 

For a fuller discussion of rational approximations for elliptic equations we 
refer the reader to [Ml 5.3.3]. 

5. RUNGE PAIRS 

The previous section is concerned with approximation on compact sets. 
We now turn to approximation on open sets. 

Let ]?i and ^2 be open subsets of the complex plane C. From Runge's 
theorem on rational approximation it follows that a necessary and sufficient 
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condition, in order that each function holomorphic in /2i can be approxi- 
mated uniformly on compact subsets of /2i by functions holomorphic in Q2-, 
is that the complement C \ J?i has no compact components in Q2- This 
'Runge' theorem was extended by Lax and Malgrange to approximation by 
solutions of elliptic equations. 

For the heat equation, which is of course not elliptic, the analogous ap- 
proximation problem was investigated by Jones [9] and Diaz [5]. 

If i? is an open set in M""*"^, we say that 12 is a Runge open set for the 
heat operator, if the entire solutions to the heat equation are dense in the 
solutions to the heat equation on Q. One has the following characterisation 
of Runge open set, proved in [9]. 

Theorem 5.1. An open set Q C M""^^ is a Runge open set for the heat 
operator if and only if, for every hyperplane H in W^^^ orthogonal to the 
time axis, H\Q has no compact components. 

Let us say that two open sets fli and fl2 in M""*"^, with i?i C /?2 5 form a 
Runge pair for the heat operator, if 5(i72) is dense in 5(i7i) in the topol- 
ogy. That is, if for each u G S{fli) there is a sequence {uj} in S{fl2)i such 
that the restrictions of the Uj to Qi and their partial derivatives of all orders 
converge to u and the corresponding partial derivatives of li, respectively, 
uniformly on compact subsets of Hi. 

The following theorem of Diaz [5] gives a necessary and sufficient condition 
that a pair (I7i, i72) of open sets in M*^"*""^, with i7i C /22) be a Runge pair 
for the heat operator. 

Theorem 5.2. Let Qi C Q2 be open sets in M""*"^. A necessary and sufficient 
condition in order that {^1,^2) be a Runge pair for the heat operator is that, 
for every hyperplane H in M""'"^ orthogonal to the time axis, the complement 
of Qi in H has no compact components in ^20 H. 

There is a small mistake in the proof of sufficiency in ^ . The assertion on 
p. 645, line 19, that an arbitrary component ol Di\Z has the form A x (a, h), 
is not correct. Perhaps this is easily fixed, however, we do not see how to 
do this. 

We denote by H{t) the hyperplane in M""*"^ which orthogonal to the time 
axis at time t. Moreover, for a set S in M""'"-'^, we denote by S{t) the slice 

snH{t). 

Example 5.3. Let Qi be the slit plane \ {(0, x) : \x\ < 1} and Q2 be the 
punctured plane M^\{(0, 0)}. Then, for t 7^ 0, the set H{t)\Qi is empty, and 
H{0) \ Qi has a single compact component, the slit {(0, x) : < 1}, which 
does not belong to the set f22 fl H{0) = {{0,x) : x / 0}. Hence it follows 
that the pair ^2) satisfies the condition of Theorem 15. 2[ However, we 
do not know whether or not ^2) is a Runge pair for the heat operator 
in M2. 

In order to highlight the problem, let us discuss some steps towards the 
proof of sufficiency in Theorem 15. 2[ They develop Theorem 3.4.3 of [7] 
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which concerns general scalar differential operators with constant coefficients 
in M""*"^. Suppose that, for each hyperplane H orthogonal to the time axis, 
every compact component of H\ni contains a component oi H\Q2- Choose 
a continuous linear functional J- on C°°(i7i) which vanishes on 5(]72)- There 
is a distribution v € £'{W^~^^) with a compact support K vn i7i, such that 
J-{u) = {v,u) for all u E C°°(/2i). Since v is orthogonal to all exponential 
solutions of the heat equation, we conclude that there exists a distribution 
g G £'{W^~^^) satisfying {—dt — A)g = v on M"+-'^. Now, the transposed 
operator —dt — Zl is hypoelliptic and g satisfies {—dt — A)g = away from 
K, so g is an infinitely differentiable function on M^^^ \ K. Moreover, the 
transposed kernel ^' satisfies ^'{—dt — A) = I on £'(W^^^), whence g = ^'*v 
on all of M"+^ Since {dt - -t',x- x')) = for ah (t, x) G Q2 

and {t' , x') G \ and all multi-indices j and a and vanishes on S{Q2), 
we get 

di,d^,g{t',x') = {v{t,x),di,d^,'P{t-t',x-x')) 
= 

for all {t', x') away from and all multi-indices j and a. Thus, g vanishes to 
infinite order on M""'"^ \ • If we prove that supp g d fix, then vanishes on 
S{Qi) and the assertion readily follows by the Hahn-Banach theorem. For 
any t, we denote by K{t) the topological hull of K{t) in Qi{t). That is, K{t) 
is the union of K{t) and all components of fii{t) \ K{t) which are relatively 
compact in {}i{t). It is always a compact subset of /2i(t). Since H{t) \ fli{t) 
has no compact components in f22{t)^ the topological hull of K{t) in fli{t) 
just amounts to the topological hull of K{t) in fl2{t)- For each t the function 
g{t, x) is an analytic function of x for x G H{t) \ K{t) (for solutions of the 
transposed heat equation are analytic in the space variables) which vanishes 
to infinite order on H{t) \ i?2(i) and for all sufficiently large x (because g has 
compact support). Denote by K the union of K{t) over all times. By the 
very construction, K \s a subset of Hi. If (i, x) K, then either x lies in the 
unbounded component of H{t) \ K{t) or x lies in the bounded component 
of H{t) \ k{t) which meets H{t) \ ni{t) and hence H{t) \ Q2{t)- In the 
former case g{t^ x) = 0, for the support of g is compact. In the latter case 
g{t,x) = 0, for g{t,x) vanishes to infinite order on H{t) \ f22{t). We thus 
conclude that the function g vanishes away from K. Unfortunately, K fails 
to be a compact subset of for it is not bounded away from the boundary 
of The arguments of [5] include approximation of g by functions with 
compact support in f?i. 

Remark 5.4. If the slices /?i(t) have no "holes," then, for each compact 
set K d Qi, the hull K is a compact subset of fli as well. Hence, the proof 
of sufficiency in Theorem \5.1\ is straightforward. 

Note that Theorem 15.21 is closely related to the Runge theorem for har- 
monic functions, for any harmonic function in a domain X C M" is a solution 
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of the heat equation in the tube domain M x <Y in , which is independent 
of i. 

6. Universality 

We shall investigate two types of universality with respect to the heat op- 
erator, universal series and universal functions. For an survey of universality, 
we refer the reader to [6j. 

Let y be a topological vector space (real or complex) and a 
sequence of vectors of V. The series 

is said to be universal in V if the set formed by the partial sums is dense in 
V. 

A subset Z" of a topological space of Baire category II is called residual 
if its complement is of Baire category I. In this topological sense, a residual 
subset of a second Baire category space contains 'most' points of this space. 

Theorem 6.1. Assume that {ujjjgN is a sequence in a metrisable topolog- 
ical vector space V over ¥ = M or ¥ = C The following statements are 
equivalent: 

1) There exists a sequence c = {cjjjgN of scalars, such that the series 
CjUj is universal in V . 

2) The set S of such sequences c is residual in endowed with the 
product topology. 

3) For each n€N, the set of finite linear combinations from {m„, m„+i, . . .} 
is dense in V. 

Proof See [H]. □ 
We remark that the space 

¥^ = l[¥, 

endowed with the product topology, is a complete metric space and hence, 
by Baire's theorem, of second category. 

Combining Theorem 16.11 with Theorem 13.21 we arrive at universal series 
in spaces of solutions to the heat equation, whose terms are simple rational 
solutions. 

Theorem 6.2. Suppose that Q is an open set in R"^^, U an arbitrary 
neighbourhood of H, and {(tj,Xj)}jgN a dense sequence in U\Q. Then there 
is a universal series for S{Q), whose terms are simple rational solutions 
with poles at {(tj, a;j)}jgp!j. That is, there exists a sequence {cj}_,GN of real 
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numbers, such that, for each u € S{Q), there is a sequence {jN}NeN, such 
that 

3n 

cj ^{t — tj,x — Xj) u 

in the topology ofC^{Q). 

Proof. According to Theorem 13. 2^ for each n G N, the set of finite hnear 
combinations of functions {^{t — tn,x — ^{t — tn+i,x — Xn+i), • • •} is 
dense in the space 5(i7i). The desired conclusion now follows directly from 
Theorem 16. li □ 



Using Examples 14.31 and 14.41 yields universal series for S{f2), whose terms 
are rational solutions with special sets of poles outside of O. 

Theorem 6.3. Suppose that f?2) is a Runge pair for the heat operator in 
]gn+i_ T/ien there is a universal series for S{fii), whose terms are elements 
of S{Q2)- More precisely, there exists a sequence {uj}j^^ in S{Q2) with the 
property that, for each solution u € S{Qi), there is a sequence {jN}N&n, 
such that 



u 



in the topology ofC°°{Qi). 



Proof. Since {fli, Q2) is a Runge pair, the subspace S{Q2) is dense in 
Since C°°{Qi) is separable and every subspace of a separable space is also 
separable, we may choose a sequence {u'j}j^^ in S[Q2), which is dense in 
S{Qi). Then, for each n G N, the sequence . . .} is also dense in 

S{Qi). By Theorem 16.11 there is a sequence {cjjjgN of real numbers, such 
that the series "^Cju'j is universal in 5(i7i). To complete the proof, we set 

Uj = Cju'j. □ 

In order to prove the existence of universal solutions, we need a lemma on 
so-called tangential approximation on an unbounded set. In the case of the 
Cauchy-Riemann operator and the Laplace operator the theory of uniform 
approximation has been developed not only on compact sets, but also on 
(possibly unbounded) closed sets. For the heat operator, we shall confine 
ourselves to very simple closed sets. Given a sequence {Kj} of disjoint 
compacta in R""*"^, we write Kj /' 00 if 

sup x)| < inf 

{t,x)eKj {t,x)€Kj+i 

inf \it,x)\ 00 
{t,x)eKj 



as J ^ 00. 
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Lemma 6.4. Let {Kj} be a sequence of convex compacta in M"^^, such 
that Kj y oo. Then, for each sequence {uj} with uj € S{Kj) and each 
sequence {ej} with Sj > 0, there is an entire solution to the heat equation 
u G 5(M"+^), such that 

sup \u{t^x) — Uj{t^x)\ < Ej 

{t,x)eKj 

for all j G N. 

Proof. For j = 1,2, . . ., choose Rj such that 

sup \{t,x)\<Rj< inf \{t,x)\. 

Set Bq = ^ and, for j = 1,2,..., let Bj be the closed ball of radius Rj 
centered at the origin in R"""*"^. We may assume that the sequence {ej} is 
decreasing. Let us proceed by induction. Since Ki is convex, by Theorem 
15.11 there exists hi G 5(M""'"^) such that \hi — ui\ < e\j2 on K\. Suppose 
functions h\,...,hj in 5(M""'"^) have been constructed with the following 
properties: 

|> hi — Uj\ < on Ki, 

(6.1) ^' 2^ 

\hi\ < ^ on Bj^i. 

Since -STj+i and Bj are convex and disjoint, by Theorem 15.11 there exists 
hj+i G 5(M"+i) such that 

l^/ii-Uj+il < ^ on Kj^i, 

i=l 

e 



Hence, for each j = 1,2,..., there exists a solution to the heat equation 
hj G 5(M''+^) which satisfies It follows that the series 

oo 

converges uniformly on compact subsets of M"^^ to an entire solution to the 
heat equation u G with the desired properties. □ 

We are now in a position to construct a universal entire solution to the 
heat equation. 

Theorem 6.5. There is an entire universal solution to the heat equation, 
that is, an entire solution u, whose translates are dense in the space of all 
entire solutions to the heat equation. Thus, for each f G 5(M"'*""^), there 
is a sequence {aj} in M""*"^, such that u{- + aj) — > /(•) in the topology of 
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Proof. Let {hj}j^^ be a sequence of entire solutions to the heat equation, 
which is dense in the space of all entire solutions to the heat equa- 

tions. Choose a sequence {Kj}j^fq of pairwise disjoint closed balls in M""*"^, 
such that Kj y oo and whose radii Rj tend to infinity. Denote by aj the 
center of Kj and set Uj{t,x) = hj{{t,x) — Oj). By Lemma [6.41 there is a 
solution u G 5(M"'~''^) such that 

sup \u{t,x) — Uj{t,x)\ < - 

for all J = 1, 2, . . .. Equivalently, 

sup \u{{t,x) + Oj) — hj{t,x)\ < - 

\{t,x)\<Rj 3 

for all j = 1,2, . . .. Since 1/j — 0, Rj oo and the sequence {/ijjjgN is 
dense in 5(M""'"^), the proof is complete. □ 

We remark that the abstract universality Theorem 16.11 shows that the 
phenomenon of universality is generic. That is, once we know the existence 
of a universal series, it turns out that 'most' series are universal. One can 
also show that the result on the existence of a universal function, Theorem 
16.51 below, is generic. 

7. Burgers' equation revisited 

We can now return to Burgers' equation (|2.2p in R""^^. We shall tacitly 
assume that the coefficient a{p) is independent of p to have an explicit 
transformation u = U{p). 

In order to be able to apply the inverse transformation p = U^^{u), we 
should remain in the domain of the inverse. It is described by 

Uq and Ui being arbitrary constants of Example 12.11 We are thus led to the 
study of solutions to the heat equation which take their values on a half- axis, 
a fixed conical set. 

Another approach we follow is to allow solutions with singularities caused 
through U. 

If / € £'{W^~^^) is an arbitrary distribution with compact support on 
M"^^, then the inhomogeneous equation 

Pt = Ap + a{p)\Vp\'^ + f{t,x), 

cf. (j2.ip . possesses a potential type solution p = U^^{<P* f) on which 
explicitly reads 

(7.1) P{t,x) = — log 1-a 



a °\ Ui 
This 'solution' no longer makes sense as a distribution on but away 

from the singularities of /. If / is the unit mass at a point (toi^o) ^ 
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we call p{t,x) and its constant multiples simple rational solutions to ()2.2p . 
More generally, by rational solutions of (j2.2p are meant solutions of the form 
()7.ip . where / is a distribution on W^^^ whose support consists of a finite 
number of points. 

From what has been proved for solutions to the heat equation we readily 
derive through the transformation u = U{p) analogous results for solutions 
to Burgers' equation. 

Theorem 7.1. For each compact set K C M"""*"^, the simple rational so- 
lutions to 12. S\) with poles outside of K are dense in the set of all smooth 
solutions on K . 

Since (j2.2p is a quasilinear equation, its solutions do not survive under 
addition. Instead, we consider compositions of solutions defined by the 
formula piop2 '■= U'^^{U {pi) + U{p2)), cf. Section [2l As usual we introduce 
the composition of infinite number of solutions, to be referred to as series, 
by 

o"LiPj = lim o . . . opiy. 

Theorem 7.2. For each open set Q in W^^^ , there is a universal series in 
the set of solutions to h2. ^) on Q, whose terms are simple rational solutions 
to 1^2. 2\) with poles outside of Q. 

That is, there is a sequence {{tj,Xj)} of points outside of ]?, and a se- 
quence {cj} in M, such that, for each solution p to (j2.2p on Q, there is a 
sequence {jAfjAfeN with 

in the topology of C°°(J?). 

Theorem 7.3. There is an entire universal solution to \2.2^) . that is, an en- 
tire solution p, whose translates are dense in the space of all entire solutions 
to / fO) . 

This universal solution is obviously given by p = U~^{u), where u is an 
entire universal solution to the heat equation whose existence is guaranteed 
by Theorem 16.51 
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